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The advanced versions of the LIGO and Virgo ground-based gravitational-wave detectors are 
expected to operate from three sites: Hanford, Livingston, and Cascina. Recent proposals have 
been made to place a fourth site in Australia or India; and there is the possibility of using the 
Large Cryogenic Gravitational Wave Telescope in Japan to further extend the network. Using 
Bayesian parameter-estimation analyses of simulated gravitational-wave signals from a range of 
coalescing-binary locations and orientations at fixed distance or signal-to-noise ratio, we study the 
improvement in parameter estimation for the proposed networks. We find that a fourth detector 
site can break degeneracies in several parameters; in particular, the localization of the source on the 
sky is improved by a factor of ~ 3-4 for an Australian site, or ~ 2.5-3.5 for an Indian site, with 
more modest improvements in distance and binary inclination estimates. This enhanced ability to 
localize sources on the sky will be crucial in any search for electromagnetic counterparts to detected 
gravitational-wave signals. 

PACS numbers: 



I. INTRODUCTION 

Gravitational waves and electromagnetic counterparts 
from the merger of compact binaries will carry comple- 
mentary information, and the successful association of 
the two types of merger signatures will allow many cru- 
cial questions in stellar and binary evolution and cosmol- 
ogy to be answered (see [U [2] and references therein). 
Good sky localization of gravitational-wave sources will 
be crucial in searching for associated electromagnetic 
transients. In this paper, we discuss ways in which relo- 
cating one of the LIGO detectors to a site in Australia 
or India could improve the prospects of multi-messenger 
gravitational-wave astronomy. 

The final S6 science run of the enhanced-LIGO gravi- 
tational wave detectors [3], along with the third science 
run of the Virgo detector [4] and GEO 600 [5], has re- 
cently concluded. Construction of the second generation 
of instruments is already underway, with 4 km Advanced 
LIGO detectors undergoing installation at the Hanford, 
WA and Livingston, LA observatories [6], with sensitiv- 
ity expected to improve by about one order of magni- 
tude. These two sites in North America are expected to 
be joined by the Advanced Virgo [7] detector, located in 
Cascina, Italy, to form a second-generation network con- 
sisting of three sites. Recently, proposals have been made 
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to add a fourth site to the network, at Gingin in West- 
ern Australia (- 31°21'28" S, 115°42'50" E) , or outside 
Bangalore in India (14°14' N, 76°26' E). The advantages 
and disadvantages of the LIGO Australia and LIGO In- 
dia proposals were studied by a working group within 
the LIGO Scientific Collaboration, and their conclusions 
for LIGO Australia reported in [S]. The Large Cryo- 
genic Gravitational Wave Telescope (LCGT) detector is a 
planned interferometer with 3 km arms and cryogenically 
cooled mirrors, located in Japan (36° 15' N, 137° 11' E) 
200 metres below ground to reduce seismic noise [5]. In 
addition to the 4-site advanced LIGO/Virgo networks, 
we also consider the possible network with two advanced 
LIGO detectors at Hanford, one at Livingston; advanced 
Virgo at Cascina and LCGT in Japan. In the following, 
we use an acronym for denoting the network based on the 
first initials of the sites involved: Australia (A), Hanford 
(H), India (I), Japan (J), Livingston (L) and Virgo (V); 
e.g., AHLV is a network consisting of LIGO detectors in 
Australia, Hanford and Livingston and the Virgo detec- 
tor. 

In this study, we make a comparison of the perfor- 
mance of the proposed networks with the performance 
of the 3-site HHLV network in terms of parameter es- 
timation for compact binaries. We focus on binary 
neutron-star systems, which are expected to be a preva- 
lent source of observable signals for the advanced detector 
network |10j . 

Using independent Bayesian analyses, we compare the 
parameter estimation performance of each network for an 
ensemble of sources spread throughout parameter space, 
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and demonstrate the improvement gained from the ad- 
dition of sites to the network. The relative improve- 
ment for individual sources is assessed at a fixed dis- 
tance or signal-to-noise ratio. Bayesian inference allows 
us to extract all of the available information about sys- 
tem parameters from the full data set. In contrast with 
the timing triangulation method discussed below, all of 
the data from all interferometers are taken into account 
coherently. And unlike Fisher information matrix tech- 
niques, which explicitly only describe the local structure 
of the parameter space and are accurate only in the limit 
of high signal-to-noise ratio (SNR), Bayesian inference 
methods search through the full parameter space, show- 
ing the global structure of posteriors, including multiple 
modes or near-degeneracies. 

Beginning with the work of several other stud- 
ies have analyzed sky-localization accuracy with different 
network configurations [l2l - fT5] . These investigations ei- 
ther used only a limited amount of information (e.g., tim- 
ing information alone) or Fisher-matrix techniques that 
may fail in a multimodal, degenerate parameter space. 
The study by Nissanke et al. [TB] used Bayesian methods 
to compare HLV, AHLV, HJLV and AHJLV networks, 
but used a single detector at Hanford for the HLV case 
(so some of their quoted improvements are simply due 
to greater SNR accessible with 4 rather than 3 detec- 
tors). Furthermore, [IB] used a population of sources 
randomly distributed in space out to z = 1, whereas we 
distributed our sources at constant distance (for nested 
sampling [NS] and Markov Chain Monte Carlo [MCMC] 
methods) or constant SNR (Fisher matrix method) and 
varied the other extrinsic parameters. This was done 
to ensure good coverage of the extrinsic parameter space 
when assessing the relative network performances. Mean- 
while, the study by Klimenko et al. focused on sky local- 
ization of transient burst sources rather than coalescing 
compact binaries with known waveforms [17] . 

The "Report of the Committee to Compare the Sci- 
entific Cases for AHLV and HHLV", which considered 
the scientific advantages of moving one of the Hanford 
detectors to Australia [5], was based in part on a pre- 
liminary version of the work presented here. Ours is the 
first study to apply Bayesian techniques to comparisons 
of parameter-estimation accuracies with a network in- 
cluding an Indian detector as well as networks with Aus- 
tralian or Japanese detectors. Furthermore, we analyze 
the impact of moving one of the detectors to Australia 
or India on the accuracy of measuring masses, distances, 
and other parameters in addition to sky localization for 
the first time [4"T] . 

This paper is organized as follows. In section II, we de- 
scribe the network configurations being analyzed. In sec- 
tion III, we introduce the analysis techniques employed 
in this paper. The simulations and their results are de- 
scribed in section IV. We conclude in section V. We also 
include an appendix discussing the impact of the change 
of network configurations on false alarm probabilities and 
detection thresholds for a fully coherent analysis. 



II. NETWORK CONFIGURATIONS 

Interferometric gravitational wave detectors are notori- 
ously bad at determining the direction of incoming radia- 
tion from short duration sources when used individually, 
as the detector has good sensitivity over a large range 
of angles. Although this allows an all-sky search to be 
performed without the need for 'pointing', it also means 
that the amplitude of the incoming signal cannot be used 
to determine its location well. To be able to resolve the 
position on the sky of a short duration gravitational wave 
source, a network of interferometers is needed. 

With more than one site it is possible to use triangu- 
lation to determine the location of an observed signal us- 
ing the observed time delay between different detectors. 
With two sites, this method can resolve the position to 
within a ring on the sky centered on the axes between 
the sites; with three sites this is reduced to two patches; 
but the addition of a fourth site allows a unique patch 
to be determined for each source. The limiting factors in 
the accuracy of the timing method are the distances be- 
tween the sites in the network and the timing accuracy of 
the sites, which is itself governed by the signal-to-noise 
ratio and the effective bandwidth of the signal in each 
detector [T^]. So, to achieve a better sky resolution for 
any particular gravitational wave, a network should con- 
sist of detectors with a wide band of sensitivity, with the 
longest possible baselines between sites. As the sensitive 
band of the noise curve is generally limited by available 
technology and fundamental noise sources, we are left 
with the option of dispersing the detectors as widely as 
possible. 

Given the necessity to build detectors on Earth, the 
maximum possible baseline would be the diameter of the 
Earth, 42.5 ms (distances converted to gravitational- wave 
travel times). The existing LIGO- Virgo network of de- 
tectors consists of three sites, at Hanford, Livingston and 
Cascina, all in the northern hemisphere. The longest 
baseline between detectors in this network, Hanford- 
Virgo, is 27 ms. With the addition of the Gingin site to 
the network, this is increased to 41 ms for the Livingston- 
Gingin baseline, close to the maximum possible, and the 
two baselines from Gingin to the other sites are all also 
of great length. The Japanese and Indian sites similarly 
increase the longest baseline, with Japan being closer to 
the U.S. sites but India closer to Virgo. All times are 
given for reference in table [i] 

In this study, we compare the networks under the as- 
sumption that all detectors are operational at the time of 
the observation of a signal. In practice, however, detec- 
tors have limited duty factors. The probability of having 
at least three non-colocated detectors up at a given time 
is higher with the larger networks than with the HHLV 
network, increasing the probability for decent source lo- 
calization. 

Although the triangulation method above captures 
well the essential reason for desiring a longer baseline, 
the methods used in this paper are based on a fully co- 
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TABLE I: Table of gravitational-wave travel times between 
sites (Hanford, Livingston, Virgo, Australia, India, Japan, 
identified by their first letters), in milliseconds. The maxi- 
mum possible baseline for a terrestrial network is 42.5 ms. 



herent Bayesian analysis of the data to extract posterior 
probability functions on the parameters of interest. This 
method naturally incorporates the information from the 
time delays between sites, but it also includes the infor- 
mation from the amplitudes and relative phases of the 
signals present in each detector. This information can be 
used to further restrict the sky location, for example by 
eliminating the secondary maximum in the sky location 
for the majority of signals in the 3-detector network. 

The addition to the network of a fourth site also gives 
a fourth separate detector orientation (instead of the 
replica of the HI detector in the HHLV network). This 
raises the possibility of improved measurement of the 
other parameters of the source, in particular, the mea- 
surement of the polarization angle %ji and inclination an- 
gle i of the gravitational wave may be expected to im- 
prove in certain cases. We perform a comprehensive com- 
parison of measurement accuracies of these and other pa- 
rameters encoded in the gravitational-wave signal from 
an inspiraling binary composed of non-spinning neutron 
stars. 



III. ANALYSES 



For this study, we used two independent Bayesian in- 
ference codes that implement two different techniques: 
the LALInferenceMCMC code (based on the SPINspi- 
ral code [T8H20] ) uses a Markov Chain Monte Carlo al- 
gorithm (MCMC) [2T], while the Inspnest code [221125] 
uses nested sampling [23] . 

Both techniques stochastically sample the parameter 
space in a search for the parameters that best match the 
observed data, simultaneously finding the set of param- 
eters that yield the best fit to the data, and determin- 
ing the accuracy of the parameter estimation. This is 
achieved by calculating the posterior probability density 
function (PDF) on the parameter space 6 of the signal, 
given the data in the frequency domain d, and a signal 



model hypothesis H, which is 
p(9\H)p(d\6,H) 



p(0\d,H) 



p(d\H) 



1 ,7 



cxp((9|#)exp ( -~{d-h{6)\d-h{6)) 



(3.1) 



where p(9\H) is the prior probability density of the pa- 
rameters 6 and h(6) is the model used to describe the 
signal in the frequency domain |25j . The noise- weighted 
residuals in the presence of Gaussian noise with power 
spectral density S n (f) are given by 



(d- h(9)\d-h(6)) = 4 



d(f)-h(f;6) 

Sn(f) 



df. (3.2) 



For these simulations, for the LIGO and Virgo detec- 
tors, including the A and I sites, we used simulated noise 
with noise power spectral density 5* n (/) similar to the 
Advanced LIGO design curve from the LIGO Algorithm 
Library (LAL) [26J. Use of the Advanced LIGO noise 
curve for the Advanced Virgo detector may change the 
absolute results slightly, however the relative improve- 
ments ought to be consistent. The noise curve fit has a 
functional form of 



Sn(f) 



111 




(3.3) 



where S = lO^Hz' 1 and f = 215 Hz. For the LCGT 
detector, an interpolated LCGT design sensitivity curve 
was used [9l [27]. 

The MCMC method used by LALInferenceMCMC ex- 
plores the parameter space with a random walk, using 
the Metropolis-Hastings algorithm to simulate samples 
from the posterior probability distribution function 
p(6\d,H). LALInferenceMCMC uses a variety of op- 
timisation techniques, including parallel tempering, to 
converge on the modes of the distribution and ensure ad- 
equate mixing of the chain [18H20| . It was started with 
randomly offset parameter values to simulate imperfect 
initial inputs from the detection pipeline, and 10 chains 
ran on each event to check for convergence. 

The nested sampling algorithm used by Inspnest oper- 
ates by generating and replacing samples from the prior 
distribution p{9), gradually shrinking the volume sam- 
pled by imposing a limit of minimum likelihood p(d\0, H) 
on each replacement sample |24j . Inspnest samples the 
prior distribution using an MCMC algorithm, which is 
optimised for exploring the structure of the limited prior, 
with jumps proposed through differential evolution, and 
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along possible degeneracies in parameter space [351 H3] • 
It is designed to compute the evidence p(d\H), but the 
output samples from the prior can be weighted appropri- 
ately and resampled to produce samples from the poste- 
rior distribution. 

The output from both Bayesian codes is a list of sam- 
ples from the posterior distribution p(0\d, h), which are 
then used to estimate the mean, variance, and percentiles 
of the distributions. We performed two-dimensional 
binning and used a greedy algorithm to compute two- 
dimensional minimum probability intervals, as explained 
in section IIVI 

As well as the two Bayesian methods, we also used the 
Fisher information matrix (FIM) to estimate the mea- 
surement uncertainties. The FIM technique, which ap- 
proximates the likelihood surface quadratically near the 
likelihood peak, has a long history in gravitational-wave 
parameter estimation (e.g., [28 ). It is known to suf- 
fer from a number of flaws, particularly in the low-SNR 
limit (when the quadratic approximation breaks down), 
or when correlations between parameters are very signif- 
icant |29j . Furthermore, the Fisher information matrix 
is entirely local, and only approximates the shape of the 
maximum where it is evaluated, ignoring other maxima 
in the global parameter space which are picked up by the 
Bayesian analyses. On the other hand, the FIM tech- 
nique is computationally inexpensive, and hence allows 
a larger number of sources to be simulated in order to 
improve statistics. 

The signal model for a non-spinning inspiral signal re- 
quires nine physical parameters: the chirp mass Ai = 
(miTO2) 3 / 5 (mi + vn-i) _1 / 5 and the symmetric mass ra- 
tio r\ = m\mi{yri\ + m2)~ 2 (where m\ and m2 are the 
individual masses), right ascension a, declination S, incli- 
nation t, orientation ip, the luminosity distance g?£, the 
time of coalescence t c and the phase at coalescence <f>. 
For the Bayesian analyses, the prior probability distri- 
bution was assumed to be isotropic on the sphere of the 
sky, and on the inclination of the binary relative to the 
line of sight, and proportional to dh' . All other priors 
are uniform unless otherwise specified. 



IV. SIMULATIONS 

For this study, we assumed that a successful detection 
has already been made, and we have the correct wave- 
form model to process the data, so we did not perform 
evidence calculations or model selection. The waveform 
model used for injection in all cases was generated by 
the LIGO Algorithm Library GeneratePPNInspiral rou- 
tine [26], which uses a time-domain approximant at 2.0 
pN order in phase and pN order in amplitude |30j . 
All injections used symmetric binary neutron star sig- 
nals with masses mi = rri2 = 1.4 as observed in 
the detector frame. For recovery and posterior calcu- 
lation, we used the frequency-domain stationary phase 
approximation TaylorF2 approximant from LAL at 2.0 



pN order (see, e.g., [3TJ[32] for additional information on 
post-Newtonian waveform approximants). 

Injections were coherently made into four network 
configurations: (i) HHLV, (ii) AHLV, (iii) HILV, (iv) 
HHJLV. All injections were performed using Gaussian 
colored noise, using the Advanced LIGO power spec- 
tral density approximated by Eq. ( |3.3[ ) for the Advanced 
LIGO and Advanced Virgo sites (including Australia and 
India), and a fit to the LCGT noise curve for the J detec- 
tor in the HHJLV network. In each of the four-site net- 
works the noise realizations were kept the same in each 
network, although the sites are moved. In the HHJLV 
network the noise realisation in the HHLV portion was 
the same as in the other 4-site networks. We also exam- 
ined the possibility of rotating the Australia detector by 
45°, and found very similar results to the AHLV case, 
which are omitted from further discussion for brevity. 

The additional A and I detectors were positioned with 
their arms oriented along the local North and East vec- 
tors projected onto the tangent plane to the Earth's sur- 
face at the site, whose locations are given in section [I] 
The J detector was aligned with its y-arm 19° from the 
local North vector as the geometry of the detector is con- 
strained by the existing underground tunnels. 

As the speeds of the methods differed, we performed a 
different set of simulations with each method, with the 
following details. 



Nested Sampling 

The nested sampling implementation was able to run 
on a large set of injections. All signals were injected at a 
luminosity distance of 180 Mpc, but the inclination an- 
gle, polarization angle, right ascension and declination 
were located ona5x6x5x5 rectangular grid in the 
a x sin 6 x cost x ip parameter space for a total of 750 
injections. This resulted in a range of SNRs between 5 
and 35, but mostly between 7 and 25. We chose low 
and high frequency cutoffs of 30 Hz and 2048 Hz, which 
included the maximum frequency of the inspiral signal 
as no merger or ring-down components were used. Of 
the 750 signals, 728 were detected in all network con- 
figurations, and it is these we will use for the summary 
statistics throughout the remainder of the paper. The 
nested sampling search, by its nature, samples the entire 
prior range given for the parameters. In this case, the 
total mass M = m\ + m 2 was assumed to be between 
2 and 35 Mq, with component masses allowed to be in 
the range mi, m2 G [1, 30] Mq. The prior distribution on 
M. used was p{M) oc d 2 L M~ 5 / 6 , chosen as an approx- 
imation to the Jeffrey's prior which sets the prior as a 
function of the Fisher matrix T on the parameter space, 

p(9) cx \J det r(0), in order to improve sampling where 
the template bank density is highest [53] , 
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MCMC 

The MCMC method was used to run on a randomly 
chosen subset of the injections that were analyzed with 
nested sampling (computational constraints prevented us 
from using the full set of injections). Results from 42 in- 
jections are used in this analysis. For consistency, the 
MCMC runs employed the same priors and frequency 
range for the overlap integral as the nested sampling runs. 
The MCMC results on individual injections matched the 
nested sampling results, therefore allowing us to gain 
extra confidence in the Bayesian parameter estimation. 
However, because of the concern that smaller numbers 
of runs could increase statistical fluctuations, we do not 
quote absolute accuracies for MCMC runs, but only com- 
pare the expected parameter estimation accuracies for 
different network configurations. 



Fisher Matrix 



F < 1, the F- width of a one-dimensional PDF was 
defined as the width of the smallest region that con- 
tains that fraction F of the posterior PDF. Thus, the 
95% width represents the width of the smallest region 
of parameter space that contains 95% of the total pos- 
terior probability. A similar approach was used for two- 
dimensional PDFs in (inclination, distance) space and 
(right ascension, declination) space, with pixels of a fixed 
size (0.25 deg 2 ) being used in a greedy algorithm to esti- 
mate the sky area for sky localization. 

We also define a "standard accuracy" (by analogy with 
the standard deviation) as the square root of the mean of 
the sum of the squared differences between the points in 
the PDF (sampled according to the posterior) and the 
true value. Thus, for a marginalized one-dimensional 
PDF, 



standard accuracy = 



\ 



1 - 
N ^ 

i=l 



(Xj 3-true) 



The FIM has the advantage of being computationally 
inexpensive, and so permits a large Monte Carlo over in- 
jections. Therefore, we used the FIM to confirm the re- 
sults of our Bayesian analyses. We used a low-frequency 
cutoff of 30 Hz and integrated up to the innermost stable 
circular orbit frequency, around 1600 Hz. 

For this study, we varied all angles in a Monte Carlo 
of 4000 points, adjusting the distance to keep the total 
network SNR equal to 30 for all injections in all network 
configurations. We chose injections with this relatively 
high value of SNR, rather than the SNR distribution used 
for nested sampling and MCMC studies, because of con- 
cerns about the the accuracy of the FIM approach outside 
the high-SNR limit. The SNRs were separately normal- 
ized to 30 for all networks, including the five-detector 
HHJLV network; therefore, results for that network from 
the FIM study are expected to be worse than those ob- 
tained with Bayesian studies, which would normally have 
higher SNRs in the five-detector network than the four- 
detector networks for a given injection. 

A. Quantities compared 



For a PDF whose mean is equal to the true value, the 
standard accuracy is just the standard deviation. For 
a delta-function PDF that is biased away from the true 
value, the standard accuracy is the error. In general, the 
standard accuracy is equal to the sum, in quadrature, 
of the standard deviation of the PDF and the difference 
between the PDF mean and the true value. 



B. Relative Improvements 

In the following sections, we present the results of each 
of the three analyses of the differences in parameter es- 
timation accuracy between the three networks. In the 
case of the two Bayesian codes, we measure the 95% con- 
fidence intervals for each parameter and compute the ra- 
tio of these in comparison to the result from the HHLV 
network for the same injection. The median value of the 
ratio is quoted in each table, along with the range incor- 
porating the 5% and 95% quantiles of the distribution. 
This gives some measure of the spread of the ratios across 
different sky positions, locations and orientations. 



Although our techniques make full nine-dimensional 
posterior PDFs available, these are unwieldy to compare 
or visualize. So we typically consider only one- or two- 
dimensional PDFs marginalized over the remaining pa- 
rameters, with examples shown in figures [4] & [3j How- 
ever, to allow us to make comparisons, we had to restrict 
ourselves to particular estimators for the PDFs, with the 
understanding that unless PDFs are extremely narrow 
or are described by a simple analytical function (e.g., a 
Gaussian) , a few estimators are not sufficient to describe 
all of the information contained in the PDFs. 

We estimated the width of a particular one- or two- 
dimensional PDF as follows. For a given fraction < 



MCMC results 

In table |Hj we average comparisons of the 95% con- 
fidence intervals across the 42 injections analyzed with 
the LALInferenceMCMC code. We show the values of 
the 95% confidence interval widths for the extended net- 
work configurations as fractions of the same widths for 
the HHLV configuration. Table [H] lists the median ratio 
of the size of the 95% probability region over all of our 
injection runs, while the 5% and 95% percentiles of the 
distribution of ratios are shown in super and subscript, 
to indicate the range of the results. 
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TABLE II: Comparative 95% interval widths and standard accuracies for one-dimensional PDFs, and comparative 95% areas 
for two-dimensional PDFs (last two lines) averaged over all injections, calculated using the MCMC algorithm. All values are 
reported as fractions of the corresponding values for the HHLV network configuration. The median values of the ratios across 
all injections are computed; the error bars correspond to the spread between the 5% and 95% quantile values of these ratios 
across all injections. See text for details. 
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Parameter 


95% width 


std. acc 


95% width 


std. acc 


95% width 


std. acc 


M 


1 


nn +0.80 


1 


00+1.92 
UU -0.56 


1.001° 


71 

47 


1 


02l 


39 

68 





Q9+0.33 
yz -0.37 





qo+0.34 
98 -0.42 


V 


1 


00 +0.78 
uu -0.38 


1 


09+1-28 
uz -0.51 


1.001° 


7(1 

44 


1 


OHo 


14 
56 





Q9+0.30 
yz -0.29 





OO+0.35 
yo -0.38 


t c 





7 o+0.54 





RO+104 
Dy -0.61 


0.69l° 


61 

46 





62l° 


90 

52 





co+0.32 
"0-0.43 





66+ ' 49 

DD -0.57 


d L 


1 


nfl +0.33 
uu -0.21 





QQ+0.15 

y8 -0.25 


1.051° 


45 
30 





911°, 


53 
30 





Q9+0.25 
9z -0.24 





qo+0.16 
98 -0.30 


a 





R7+0.58 
"' -0.48 





ci +0.88 
D1 -0.55 


0.62l° 


52 
47 





56l 


00 
50 





60+ 040 

""-0.44 





c: fi +0.59 
°"-0.52 


8 





cn+0.70 
ou -0.40 





OQ+1.12 
oy -0.34 


0.611° 


59 
48 





48l 


02 
41 





C7+0.33 
"' -0.49 





c:q+0.53 
,Jy -0.49 


1 





nq+0.24 
3 °-0.46 





Q1 +0.35 
9i -0.58 


0.86l° 


2* 
47 





83l° 


56 
56 





Qc+0.19 
8"_0.44 





OD+0.26 
9u -0.59 


4, 


1 


nn+0.19 
uu -0.52 





qq+o.32 

yy -0.51 


0.97l° 


11 
66 





9H°o 


32 
63 





OO+0.11 
yo -0.52 





O7+0.22 
y ' -0.56 


a — 8 





O.9+0.78 
JZ -0.24 






0.431°, 


87 
32 











OO+0.51 
°°-0.22 






d L — l 





qn+0.29 
au -0.45 






0.88l° 


44 
49 











ro+0.24 
do -0.27 







TABLE III: Comparative 95% probability interval widths and standard accuracies between alternative network configurations 
and the HHLV network, calculated using the nested sampling algorithm. The median values of the ratios across 728 injections 
detected in all networks are computed. Errors quoted correspond to the 5% and 95% quantiles of the distribution of ratios, as 
in table [Tl] For a — 8 plane, 95% probability intervals are calculated using a greedy binning algorithm. 



Parameter 


AHLV/HHLV 


HILV/HHLV 


HHJLV/HHLV 


M 


1.00 


1.00 


1.02 


V 


1.00 


1.00 


1.02 


t c (msec) 


0.54 


0.56 


0.65 


<Il 


0.77 


0.66 


0.71 


a (deg) 


0.66 


0.56 


0.50 


5 (deg) 


0.31 


0.57 


0.62 


t (deg) 


0.76 


0.67 


0.70 


tp (deg) 


0.77 


0.67 


0.69 


<t>o (deg) 


0.88 


0.84 


0.91 


a — 8 (deg 2 ) 


0.29 


0.35 


0.47 



TABLE IV: Ratios of median standard deviations for each parameter and the sky area as reported by the Fisher information 
matrix, as a function of network configuration. All injections are re-normalized to an SNR of 30 in the given network, so the 
last column, corresponding to a five-detector network re-normalized to the same SNR as the four-detector networks, can not 
be directly compared to similar columns in the tables above. 

As discussed in more detail below, the accuracy with cantly affected by the network configuration. This con- 
which individual masses can be recovered is not signifi- forms to our expectation that mass measurements, which 
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FIG. 1: Cumulative histogram of the 95% confidence interval 
for the area of the sky in square degrees, estimated using the 
Inspnest analysis of 728 signals which were detected in all net- 
work configurations, covering the range of sky locations and 
orientations of the binary. Different lines correspond to differ- 
ent networks: green dot-dashed line, AHLV; red dashed line, 
HHLV; cyan solid line, HILV; black solid line HHJLV. The ex- 
tended networks show significantly better performance than 
the 3-site network, with AHLV offering the highest fraction of 
signals resolved to better than 10 square degrees. HILV and 
HHJLV show similar performance to each other for a large 
fraction of the signals, but the HHJLV network avoids the 
tail of poorly resolved signals located in regions of parameter 
space with poor sensitivity in HILV. 



come from waveform phasing, are constrained primar- 
ily by the total network SNR. Sky localization accuracy 
can be significantly improved in both directions when a 
fourth site is added to the network. Timing accuracy at 
the geocenter is strongly correlated with sky localization 
and is similarly improved. A fourth site also moderately 
improves the accuracy of distance measurements. 

We should particularly point out the next-to-last line 
of the table, "a — 5". The area of this 2-dimensional 
PDF is a direct measure of the uncertainty in estimat- 
ing the position of the source on the sky. The error box 
shrinks by a factor of ~ 4 when the second Hanford detec- 
tor is moved to Australia or India because of the much- 
improved north-south baseline. This improvement in sky 
localization accuracy will make the detection of an elec- 
tromagnetic counterpart to the gravitational-wave source 
more feasible, and is perhaps the biggest boon in moving 
one of the Hanford detectors to India or Australia. 



Inspnest results 

In figure [l] we show the results of determining the sky 
location of 728 injections with the Inspnest nested sam- 
pling code. The lines show the cumulative histograms 
for the 95% confidence interval on the sky location of 
the binary for each of the network configurations consid- 
ered. As expected, the addition of a fourth site to the 
network yields a significant improvement in the resolu- 
tion of the sky location, due to the increased baseline 



for relative timing measurements. The greatest improve- 
ment is seen for the addition of the the Australian de- 
tector, yielding the longest north-south baseline. HILV 
and HHJLV are both great improvements on the HHLV 
network, and have similar performance for the majority 
of signals, although HHJLV suffers from fewer outliers 
with sky resolution greater than 100 square degrees, pos- 
sibly due to the increased signal-to-noise ratio from the 
additional detector. 

This improved performance is reflected in the compari- 
son of the median 95% confidence interval for each of the 
networks, given in table [V] The med ian resolution for all 
sources shrinks from 30.25 square degrees to 6.625 square 
degrees between HHLV and the most accurate AHLV net- 
work. In addition to the reduction of the median resolu- 
tion, these improvements are also reflected in the relative 
resolutions of the networks. Table IIIII shows the median 
ratio of the 95% probability region size for the extended 
networks against HHLV; note that the median of the ra- 
tios for the sky location intervals (a — S row) is not the 
same as the ratio of the medians in table [Vj 

The qualitative behavior inferred from the LALInfer- 
enceMCMC and Inspnest runs is in excellent agreement. 
Some differences between the "typical" results tables |III| 
and [TT] are due to the impact of statistical fluctuations in 
the smaller subset of injections analyzed with LALInfer- 
enceMCMC ; however, we have confirmed quantitative 
agreement between the two methods on the individual 
injections being analyzed. 

For a more detailed look at the performance of the 
networks for the non-sky-location parameters, we show 
in figure [5] the cumulative probability distributions for 
the fraction of signals found within a given probability 
interval width for M, r], <1l, l, ip and t c . Although our 
range of injections at constant distance is astrophysically 
unrealistic, these figures give a good idea of the relative 
performance of the networks across an isotropic distribu- 
tion of sky position, polarisation and inclination angles. 
We can see immediately that the improvement in the res- 
olution of the chirp mass and r\ parameters is marginal, 



as evident from tables III The slightly improved perfor- 
mance for HHJLV can be attributed to the higher signal- 
to-noise ratio in that network thanks to the fifth detector. 

With the distance and inclination angle parameters, 
we see slight improvements, with the biggest effect again 
being produced by the additional SNR in the HHJLV 
network. Due to degeneracies in the parameter space, 
the effect is more pronounced when the two dimensional 
distribution is considered, as in section [IV C| From table 
III[ we see that the median relative improvement of the 
polarization angle ip resolution is minimal, however in fig- 
ure [5] it is apparent that for the 4-site networks there are 
a greater fraction of injections which are resolved well. 
This is explained by the large fraction of signals (~ 50%) 
which have a very broad distribution where the ip pa- 
rameter is degenerate with the phase of the signal, when 
l is close to or ir. Although the median of the im- 
provement ratio for individual sources is ~ 1, the mean 



is - 0.95 for HALV and - 0.85 for HILV and HHJLV, 
indicating that the sources which can be resolved in ip are 
better resolved in the 4-site networks. The t c distribu- 
tion shows similar performance for all 4-site networks for 
most signals, with slightly fewer signals poorly resolved 
in time with the HILV network than AHLV, reflecting 
the shorter tails of the distributions in figures [T] and [2] 
(sky localization is strongly correlated with the accuracy 
of timing the gravitational wave at the geocenter). 



Network 


Median 95% conf. int. 


AHLV 
HHLV 
HILV 
HHJLV 


6.625 deg 2 
30.25 deg 2 
9 deg 2 
9.5 deg 2 



TABLE V: Median 95% confidence intervals in square degrees 
for each network configuration. 



Fisher matrix results 

The Fisher information matrix results are presented in 
table HVl These are ratios of the median standard de- 
viations for each parameter, taken from a Monte Carlo 
over angles (with distance adjusted to keep the total net- 
work SNR equal to 30). According to the FIM results, 
the greatest change in sky localization comes from the 
improvement in the measurement of the declination an- 
gle with a network that includes an Australian detector. 
The declination angle can be measured more accurately 
by a factor of ~ 3 with an AHLV network relative to the 
HHLV network, as a detector in the Southern hemisphere 
greatly improves the latitudinal baseline of the network, 
allowing for superior angular resolution in that direction. 
The last line in tabic IIVI indicates the Fisher-matrix es- 
timate of the sky area, defined following Eq. (43) of [31] . 
Figure [2] shows a histogram of the sky-area accuracy for 
two network configurations. 

C. Parameter Degeneracies 

In cases where a signal is observed with marginal SNR 
in one or more detectors, or when the orientation of the 
binary is unfavourable (i.e. l is close to 90°), a degen- 
eracy between two sky locations often emerges. The ad- 
dition of the fourth site to the network ensures that the 
source can be localised to a single region on the sky, as 
shown in the 2-dimensional PDFs of figure [4] which com- 
pares the 95% probability regions for different networks 
but the same injection. Note the two areas on the sky 
in the left panel, indicating a degeneracy that exists in 
a network with three sites but is removed when a fourth 
site is added. 

Moderate improvements also appear in inclination 
and luminosity distance measurements. Here, the one- 
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FIG. 2: A cumulative histogram over injections of the esti- 
mated accuracies in sky localization obtained via the Fisher 
information matrix technique. Blue, green, red and cyan 
are the distributions of sky-localization uncertainties for the 
AHLV, HHLV, HHJLV and HILV networks, respectively. All 
injections are independently normalized to an SNR of 30 for 
each network configuration, to avoid concerns about the trust- 
worthiness of FIM results at low SNRs. The relative shapes 
of the histograms are thus more relevant than their actual 
values. 



dimensional estimators used in Tables [TTJ |III| and |IV| do 
not tell the full story because of the strong correlation 
between inclination and luminosity distance. The sec- 
ondary maximum on the sky position also corresponds to 
a secondary mode in the cos i parameter corresponding 
to the transformation l t— > n — i. In fact, the addition of 
a fourth site allows this degeneracy to be broken, as can 
be seen by comparing the marginalized one-dimensional 
PDFs for a sample injection the HHLV and AHLV net- 
work configurations, plotted in figure [3| Although the 
width did not shrink significantly from red to blue pos- 
teriors, with the AHLV network configuration the poste- 
riors are unimodal and centered on the true values, as a 
degeneracy in the inclination-distance space is broken. 

Both these effects are clearly visible in figure |4j which 
shows the breaking of the degeneracy in d^-i-a-d space 
as projected onto the d^-i and a— 5 planes. In this exam- 
ple all the expanded networks allow the secondary max- 
imum to be eliminated. Figure [6] shows the cumulative 
distribution of the 95% probability region for the d^-t 
space, with the improvement in two dimensions more 
prominent than looking at d^ or i individually (c.f. figure 

On the other hand, the accuracy with which mass 
parameters are measured does not improve as a conse- 
quence of moving the site of the fourth detector to India 
or Australia. We can speculate that the reason for this is 
that information about masses comes primarily from the 
phase evolution of the signal, and the accuracy in mass 
estimation is predominantly set by the overall network 
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FIG. 3: Comparison of the one-dimensional PDFs for a typi- 
cal source as detected by the HHLV network (red) and AHLV 
network (blue). Note the bimodal posteriors in right ascen- 
sion and declination for the HHLV network vs. unimodal ones 
for the AHLV network. The latter network also allows for bet- 
ter estimates of the posteriors for inclination and luminosity 
distance, which is not properly reflected by the simple estima- 
tors of the PDF width used in table [H] Dashed lines indicate 
the true injected values (different true values of the luminos- 
ity distance were used for the HHLV and AHLV injections so 
that the total network SNR is 15 in both cases). 



SNR. Meanwhile, masses do not strongly correlate with 
extrinsic parameters (with the exception of the time of 
coalescence), so their estimation is not significantly im- 
proved by better sky localization or inclination measure- 
ments. 

The results from the Fisher information matrix are in 
qualitative agreement with the two Bayesian approaches 
regarding the partial breaking of the distance/inclination 
degeneracy achieved by moving a detector to Australia or 
India (leading to marginal improvements in both param- 
eters, see table IV). They also indicate that the accuracy 



with which masses can be measured is not affected by 
the network choice. 



V. CONCLUSIONS 

In this paper we studied the effect on parameter es- 
timation of different networks of advanced detectors. 
We employed two different Bayesian techniques and the 
Fisher information matrix to estimate the accuracy of 
parameter recovery. We analysed a set of injections dis- 
tributed in a grid in the extrinsic parameter space (with- 
out varying the mass and distance of injections) with 
the Inspnest code, and verified the results with LAL- 
InferenceMCMC . We performed a large scale Monte 
Carlo simulation using the Fisher matrix method with 
constant-SNR injections. We found consistent results be- 



tween the three methods, pointing to significant gains in 
sky localization (typically by a factor of ~ 3 — 4) and 
modest gains in distance and inclination measurements 
with a network including a fourth site. We found that the 
4-site networks are able to better resolve the polarisation 
angle of the source, in the cases where this is possible. 
We found no significant effect on mass measurements. 

Comparing the different network configurations, we 
found, as expected, the strongest improvement in sky lo- 
calization capability when the longest baseline (namely 
AHLV) was used, but that a site in India also provides a 
significant improvement in sky resolution. The HHJLV 
network, with the shortest extra baseline, provides the 
weakest improvement in sky resolution at a fixed signal- 
to-noise ratio; however, the fifth detector in this network 
can mitigate this, for an overall performance similar to 
HILV, but with fewer signals in the tail of the distribution 
with poor resolution. 

We also find good agreement with previous work. In 
particular, Fairhurst [T3] finds 20-50% of signals localised 
within 20deg 2 for HHLV, and up to 20% within 5deg 2 
with HHJLV, for an ensemble of sources at fixed distance 
of 160 Mpc, in good agreement with figure [TJ Despite the 
use of a different population of sources, Nissanke et al 
[TB] find results which qualitatively agree with our own. 
Comparison of the Fisher matrix results in figure [2] with 
Wen and Chen [T5] shows good qualitative agreement 
with their expected distribution for the HLV network at 
fixed SNR of 15, when taking into account a factor of 
(30/15) 2 = 4 for the difference in SNR used (30 in our 
case, 15 in theirs). 

In the present study, we focused on binary neutron 
stars (NS), which are the most "confident" source for the 
advanced detectors, but which are not expected to have 
significant spins [35]. On the other hand, black holes 
(BH) in NS-BH or BH-BH binaries can be rapidly spin- 
ning. Previous studies (see, e.g., [TS]) have shown that 
the presence of spin in one or both binary components 
can aid sky localization by providing additional polar- 
ization information through the precession inherent in 
misaligned spinning binaries. Localization may be fur- 
ther enhanced when a signal from a spinning binary is 
captured by a four-detector network; on the other hand, 
improved resolution of extrinsic parameters with the help 
of a fourth detector site may aid in the reconstruction of 
astrophysically interesting quantities such as spin-orbit 
misalignment angles. 

The improved ability to localize sources on the sky will 
be crucial in any search for electromagnetic counterparts 
to detected gravitational- wave signals (e.g., [TJ El [36]). 
Accurate measurements of the location of the merging 
binary can also be useful even in the absence of electro- 
magnetic counterparts, for example, in measuring typ- 
ical binary kick velocities [37]. We thus conclude that 
scientific considerations strongly favor an international 
gravitational wave network with four or more sites. 
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FIG. 4: An example of the ability of a four-site network to break the degeneracy between multiple modes in a — 8 and cos i—dt 
parameter space. This example shows the 95% confidence intervals from the HHLV (purple), AHLV (blue), HILV (green) 
and HHJLV (orange) networks. The sky position can be confined to one region with the four-site network, while the partial 
inclination angle degeneracy upon reflection is broken. The location of the injection in parameter space is indicated with a 
star. 
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Appendix: False alarm probabilities in a coherent 
search 



It is instructive to ask how different network configu- 
rations affect false alarm probabilities (FAPs), or, alter- 
natively, what the detection thresholds would need to be 
for a fixed FAP. We consider only the case of a fully co- 
herent search; the FAP of triggers in searches that only 
require coincidence between detectors in mass and time 
of coalescence parameters, such as [38l[39], should not de- 
pend on the accuracy of measuring the sky location and 
other extrinsic parameters. If a coherent search is used, 



and all detectors are assumed to have independent noise, 
we might expect that the coherence constraint would be 
stronger for two co-located detectors that should have 
the same signal, and therefore HHLV should have a lower 
FAP than AHLV or HILV for a fixed network SNR. 

We provide a Bayesian treatment of this question by 
comparing the odds ratio between the coherent signal 
hypothesis GW and the noise hypothesis N for different 
network configurations. The odds ratio is just the ratio 
of posterior probabilities for the two models, 



D, 



GW,N 



P(GW\d) _ P{GW) P(d\GW) 
P(N\d) ~ P{N) P(d\N) ' 



(A.l) 



The only term that depends on the network configura- 
tion is the evidence for the presence of a signal, Zqw = 
P(d\GW), which can be written as 



Z GW = / P (0\GW)p(d\0,GW)dO. 



(A.2) 



We assume for the sake of simplicity that the prior 
p(6\GW) — k is constant in the small region where the 
likelihood is significant, and the likelihood p(d\0,GW) 
is Gaussian in the model parameters about a maximum 
L max at 9q (as assumed for the Fisher matrix calcula- 



tion) , 



p(d\9,GW) « L max exp 



(A.3) 
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This yields 

Z G w = fcL max (2 7 r)- Ar / 2 VdctC, (A.4) 

where C is the covariance matrix of the parameters, N 
is the dimensionality of the model, and L max is the max- 
imum likelihood of the data, given by 

L max cx exp (d-~h{6 Q )\d- h(e Q )^j = ^. (A.5) 

For a fixed network SNR (maximum likelihood), the ev- 
idence is therefore proportional only to \/det C, which 
scales with the size of the region in parameter space to 
which a signal's parameters can be constrained. 

Thus, for a fixed SNR, the Bayes factor is larger when 
parameters are less precisely estimated. Due to the 
shrinking size of the posterior distribution in sky loca- 
tion, distance and inclination, we might expect a decrease 
by a factor of ~ 5 in the allowed fraction of the prior vol- 
ume between the HHLV and AHLV or HILV networks. 
Then, for a fixed detection threshold p m i n , the odds ra- 
tio will be a factor ~ 5 smaller for the AHLV or HILV 



network than for the HHLV network. Conversely, if we 
want to keep the same false alarm probability (same mini- 
mal odds ratio required for detection) , the SNR threshold 
p%un ~ 21ogL max for the AHLV or HILV network must 
increase by 2 log 5 relative to the HHLV network. This 
corresponds to a very modest increase in the SNR: for 
example, if the threshold is 12 for the HHLV network, it 
will only rise to 12.13 for the AHLV or HILV network. 

The analysis above assumes that the noise in all in- 
terferometers is uncorrelated. While that is a reasonable 
assumption for distant interferometers, various environ- 
mental factors can lead to correlated noise in co-located 
interferometers, such as the two possible detectors at 
Hanford. This correlation can increase the threshold nec- 
essary for detection. For example, at times when only 
two of the four detectors in our presumed network are 
operational, it may not be possible to make a convincing 
detection if only the two Hanford detectors are operating 
from the HHLV network, while any two detectors from 
the AHLV or HILV network have a chance to detect a 
sufficiently loud signal. 



[1] J. S. Bloom et al. (2009), 0902.1527. 

[2] L. K. Nuttall and P. J. Sutton, Phys. Rev. D 82, 102002 

(2010) . 

[3] B. P. Abbott et al. (The LIGO Scientific Collaboration), 
Reports on Progress in Physics 72, 76901 (2 009), URL 
|http : //stacks . iop . org/0034-4885/72/i=7/a=07 6901 

[4] F. Acernese et al. (The Virgo Collaboration), Classical 
and Quantum Gravity 25, 184001 (2008), URL |http : | 
//stacks . iop . org/0264-9381/25/i=18/a=184001| 

[5] H. Grote and the LIGO Scientific Collaboration, Classi- 
cal and Quantum Gravity 25, 114043 (2008). 

[6] G. M. Harry and the LIGO Scientific Collaboration, Clas- 
sical and Quantum Gravity 27, 084006 (2010). 

[7] Virgo Technical Report VIR-0027A-09, Virgo (2009). 

[8] S. Finn, P. Fritschel, S. Klimenko, F. Raab, 
B. Sathyaprakash, P. Saulson, and R. Weiss, LIGO 
Document LIGO-T1000251, The LIGO Scien- 
tific Collaboration (2010), https://dcc.ligo.org/cgi- 
bin/DocDB/ShowDocument ?docid= 1 1 604 . 

[9] K. Kuroda (The LCGT Collaboration), Class. Quantum 
Grav. 27, 4004 (2010). 
[10] J. Abadie et al. (The LIGO Scientific Collaboration), 

Classical and Quantum Gravity 27, 3001 (2010). 
[11] P. Jaranowski, K. D. Kokkotas, A. Krolak, and 
G. Tsegas, Classical and Quantum Gravity 13, 1279 
(1996). 

[12] S. Fairhurst, New Journal of Physics 11, 123006 (2009). 
[13] S. Fairhurst, Classical and Quantum Gravity 28, 105021 

(2011) . 

[14] B. F. Schutz (2011), 1102.5421. 

[15] L. Wen and Y. Chen, Phys. Rev. D 81 (2010), 1003.2504. 
[16] S. M. Nissanke, J. L. Sievers, N. Dalai, and D. E. Holz 

(2011), 1105.3184. 
[17] S. Klimenko, G. Vedovato, M. Drago, G. Mazzolo, 

G. Mitselmakher, C. Pankow, G. Prodi, V. Re, F. Salemi, 



and I. Yakushin, Phys. Rev. D 83, 102001 (2011), 
1101.5408. 

[18] M. van der Sluys, V. Raymond, I. Mandel, C. Rover, 
N. Christensen, V. Kalogera, R. Meyer, and A. Vec- 
chio, Classical and Quantum Gravity 25, 184011 (2008), 
0805.1689. 

[19] M. V. van der Sluys, C. Rover, A. Stroeer, V. Ray- 
mond, I. Mandel, N. Christensen, V. Kalogera, R. Meyer, 
and A. Vecchio, Astrophysical Journal 688, L61 (2008), 
0710.1897. 

[20] V. Raymond, M. V. van der Sluys, I. Mandel, 
V. Kalogera, C. Rover, and N. Christensen, Classical and 
Quantum Gravity 26, 114007 (2009), 0812.4302. 

[21] W. R. Gilks, S. Richardson, and D. J. Spiegelhalter, 
Markov chain Monte Carlo in practice (Chapman & 
Hall/CRC, 1996). 

[22] J. Veitch and A. Vecchio, Classical and Quantum Gravity 
25, 184010 (2008). 

[23] J. Veitch and A. Vecchio, Phys. Rev. D 81, 062003 
(2010). 

[24] J. Skilling, Bayesian Analysis 1, 833 (2006). 

[25] C. Rover, R. Meyer, and N. Christensen, Classical and 

Quantum Gravity 23, 4895 (2006), arXiv:gr-qc/0602067. 
[26] LSC Algorithm Library software, 

http : //www. lsc-group .phys .uwm . edu/lal. 
[27] K. Kuroda (2011), Private communication. 
[28] C. Cutler and E. Flanagan, Phys. Rev. D 49, 2658 (1994). 
[29] M. Vallisneri, Phys. Rev. D 77 (2008). 
[30] L. Blanchet, B. R. Iyer, C. M. Will, and A. G. Wiseman, 

Class. Quantum Grav. 13, 575 (1996). 
[31] A. Buonanno, B. R. Iyer, E. Ochsner, Y. Pan, and 

B. S. Sathyaprakash, Phys. Rev. D 80, 084043 (2009), 

0907.0700. 

[32] T. Damour, B. R. Iyer, and B. S. Sathyaprakash, 
Phys. Rev. D 63, 044023 (2001). 



12 



[33] H. Jeffreys, Proc. R. Soc. Lond. A 186, 453 (1946). 
[34] L. Barack and C. Cutler, Phys. Rev. D 69, 082005 (2004). 
[35] I. Mandel and R. O'Shaughnessy, Classical and Quantum 

Gravity 27, 114007 (2010), 0912.1074. 
[36] I. Mandel, L. Z. Kelley, and E. Ramirez-Ruiz, ArXiv e- 

prints (2011), 1111.0005. 
[37] L. Z. Kelley, E. Ramirez-Ruiz, M. Zemp, J. Diemand, 

and I. Mandel, Astro-physical Journal 725, L91 (2010), 

1011.1256. 

[38] J. Abadie et al. (LIGO Scientific Collaboration and 
Virgo Collaboration), Phys. Rev. D 82, 102001 (2010), 
1005.4655. 



[39] J. Abadie et al. (LIGO Scientific Collaboration and 

Virgo Collaboration), Submitted to PRD (2011), 

arXiv:1102.3781. 
[40] S. Nissanke, D. E. Holz, S. A. Hughes, N. Dalai, and 

J. L. Sievers, Astrophysical Journal 725, 496 (2010), 

0904.1017. 

[41] The impact of an Australian detector on distance and 
inclination measurements was discussed in [40] in a very 
different context, where the sky location was assumed to 
be known perfectly because of a presumed detection of 
an electromagnetic counterpart. 



13 




FIG. 5: Cumulative histograms comparing the ability of the networks to resolve the individual parameters of main interest. 
Plots show fraction of events found within a given width of the confidence interval. The chirp mass M and symmetric mass 
ratio r; (top row) do not benefit significantly from a 4th site, except through better uniformity of SNR across the sky. The 
distance and inclination angle (middle row) show slight improvements in their resolvability but the effect is most marked in the 
2D confidence intervals (see figure RjJ. The polarization angle ip shows interesting behaviour, with the additional sites able to 
resolve a greater fraction of the signal to within a given interval of tp, but with the median ratio of intervals being approximately 
unity (see table III I. This is explained by the sizeable fraction of signals where the polarization angle is not well resolved. The 
resolution of the time of coalescence parameter t c is similar for all 4-site networks, with 90% of signals resolved within 5 ms. 
The longer tail of less well resolved times agrees with figures [l] and [2] which show a hint of longer tail in sky resolution for 
AHLV compared to HILV. 
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FIG. 6: Cumuative histogram comparing the ability of the 
networks to resolve the signal in 2D cLl—i space. The im- 
provement in parameter estimation with the 4-site networks 
due to elimination of the secondary mode is more striking 
here than in the individual and l parameters (c.f. figure 



